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We investigate the structure of the instantaneous Bethe- 
Salpeter equation for gq-bound states in the general case of 
unequal quark masses and develop a numerical scheme for the 
calculation of mass spectra and Bethe-Salpeter amplitudes. In 
order to appreciate the merits of the various competing mod- 
els beyond the reproduction of the mass spectra we present 
explicit formulas to calculate electroweak decays. The results 
for an explicit quark model will be compared to experimental 
data in a subsequent paper Q. 



I. INTRODUCTION 

Despite many efforts the bound state problem of QCD 
is still far from being well understood. One of the main 
tasks is to investigate the relevance of quarks as dynam- 
ical degrees of freedom in hadronic bound states. Since 
a relativistic treatment for the quarks in deeply bound 
states is essential, the Bethe-Salpeter(BS)-equation [^|^ 
provides a suitable starting point. Unlike in QED one 
cannot use perturbation theory to obtain useful approxi- 
mations for the interaction kernel in QCD. Therefore our 
knowledge of the interaction between quarks is still quite 
fragmentary and various phenomenological alternatives 
have to be tested. In the present work we will restrict 
ourselves to gq-states. The use of general gg-interaction 
kernels depending on the relative time variable leads to 
serious conceptional and practical problems [^j5| . There- 
fore it is very useful at this point to make the simplify- 
ing assumption that the BS-kernel can be approximated 
by an effective interaction that is instantaneous in the 
rest frame of the bound state. The BS-equation then re- 
duces to the (full) Salpeter equation ^ that has been 
investigated for gg-states e.g. by Llewellyn Smith M, Le 
Yaouanc and coworkers |^ and recently by Lagae . 

In the present paper we analyze the properties of the 
Salpeter equation for thegeneral case of unequal quark 
masses m Sec.||. In Secg we present a flexible numer- 
ical treatment of this equation based on the variational 
principle of ref. |9). We use the block structure of the 
Salpeter amplitude to derive an expansion in terms of 
a complete set of basis functions which leads to a ma- 
trix equation analogous to the RPA (random phase ap- 
proximation) equation well known in nuclear theory pO[ . 
Our numerical method can be applied to a wide class 
of phenomenological interaction kernels. In Sec. IV we 



show how to reconstruct the full BS-amplitudes from the 
Salpeter amplitudes and present the calculation of some 
decay observables within the Mandelstam formalism in- 
cluding the decay tt'^ — > 27. Concluding remarks given 
are m SecJ^. 

The method obtained in this paper will be applied to 
an explicit quark model for light mesons in a subsequent 
paper M. 



II. GENERAL PROPERTIES OF THE SALPETER 
EQUATION 

A. Formulating the Salpeter equation 

The BS-amplitude x for a fermion-antifermion bound 
state \P) is defined by 

[XpixuX2)]^^^{0\T^„{x,)^^{x2)\ P) (1) 

where P is the four-momentum of the bound state, T 
denotes the time ordering for the fermion operators ^ 
and a, /3 stand for spinor, flavor and color indices. Due to 
translational invariance the Fourier transformation can 
be written as 



XP{XI,X2) 



-iPX 



(2^ 



''xp{p) 



(2) 



where xi = X + ri2X, X2 = X — rjix with the conjugate 
momenta pi — rjiP + p, p2 — r]2P — p. Here rji, 772 are 
two arbitrary real numbers satisfying 771-1-772 = 1- The 
BS-equation for xp{p) then reads 



xp{p) 



(3) 



S({Pi) 



[-iK{P,p,p')Xp{p')]S^{-p2) 



also shown in graphical form in Fig.^. The interaction 
kernel K{P,p,p') generally acts on xp{p') as 

[K{P,p,p')xp{p')]af3 = ^ K{P,p,p')aa\l3l3'Xp{p')a' (}' 
a' P' 

(4) 

For an interaction that is instantaneous in the rest frame 
of the bound state with momentum P = (Af , 0) the BS- 
kernel can be written as 



K{P,p,p') 



P={M,0) 



V{p,p') 



(5) 
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which can also be formulated in a covariant way as 



(6) 



where p± = p — [Ppj P^) P is perpendicular to P. In 
practical calculations one has to justify this ansatz a pos- 
teriori by investigating its consequences in the framework 
of explicit models. 

Furthermore we will approximate the full quark 
propagators S^{p) by bare propagators Sf{pi) ~ 
i {jPi — nrii + ie) ^ where nii and m2 are interpreted as 
effective masses for the quark and the antiquark. This ap- 
proximation has been criticized |]lT| because free propa- 
gators might be incompatible with a confining kernel. On 
the other hand one can argue that this choice naturally 
leads to nonrelativistic potential models that have been 
applied successfully to heav y q uarkonia (a recent model 
calculation is presented in ||12| , for an extensive review 
see [|l3|). We thus feel that free propagators should be a 
reasonable effective parameterization at least for heavy 
quarks. It is still an open question whether free propa- 
gators can also be applied to light quarks, and one has 
to investigate this problem within explicit models. 

With an instantaneous BS-kernel and bare propagators 
with effective quark masses one can perform the p'^ inte- 
grals in the BS-equation in the rest frame of the bound 
state with mass M (see e.g. [Q) and thus arrives at the 
(full) Salpeter equation 



Hp) 



dy A^ip) 70 [{Vip,p') Hp')] 7° A+i-p) 



(27r)3 M + LLI1+ UJ2 

d^p' A+ (p) 70 [{y{p, p') Hp')] 7° A- i-p) 



(27r)3 



AI — LUl — UJ2 



(7) 



with LOi = \/p^ + rnf and the projection operators 
Af{p) = {uji ± Hi{p))/{2uji) on positive and negative en- 
ergies. Here Hi{p) — j'^{'yp + mi) is the standard Dirac 
hamiltonian. We also have introduced the Salpeter am- 
plitude $ by 



Hp) 



dp"xp{p',p) 



P=(ALO) 



(8) 



For weakly bound states with \p\/mi ^ 1 and M 
TOi -f TO2 one has 



1 



1 



(9) 



so that the first term in eq.(0) can be dropped. This 
leads to the so called reduced Salpeter equation, which 
has been used in various studies of relativistic bound 
states (see e.g. the work of Gara et al. and refer- 
ences therein). In the case of light quarks, however, the 
use of the reduced Salpeter equation is dubious, espe- 
cially for deeply bound states like the pion. Quark mod- 
els for light quarks should therefore be based on the full 
Salpeter equation eq.(0). 



Let a, /9 in eq.(^ refer to Dirac indices in the standard 
Dirac representation of ref. Then $ is a 4x4-matrix 
in spinor space that can be written in block matrix form 
as 



$ = 



(10) 



where each component is a 2x2-matrix. Applying A^{p) 
from the left hand side and AJ(— p) from the right hand 
side to the Salpeter equation leads to 



At[p)Hp)At{^p) = Q 
A^{p)Hp)A^i-p) = 



(11) 



These relations allow us to express 3>+ , 3> + in terms 
of $++, $ as 



= +Cl<I>++S - C2S$ 

= -ci<I> s + C2SH 



(12) 



with the shorthand notation s — ap, Ci = U}i/{ujim2 + 
102^1). We thus find that $ can be written as 



$ $($^ 



(13) 



with $ being a bilinear function. One can interpret $++ 
as the upper component and as the lower component 
of $, as can be seen in the nonrelativistic limit where 
vanishes for solutions that fulfill M « mi-\-m2 and where 
ia2 becomes the usual Schrodinger wave function. 
For further discussion it is useful to rewrite the 
Salpeter equation in the form of an eigenvalue problem 
for the bound state mass M . We follow the treatment of 
Lagae Q and define 

^{p) :=<i>(p)7° 
[W{p,p') V'(p')] := 7° [V{p,p') Hp')] 
The Salpeter equation can now be written as 
(HV^)(p) = MV(p) 

where 

{n^){p) = H,{p)i,{p) - mH2{p) 
d^p' 



(14) 
(15) 

(16) 



(2^)3 
(2^) 



A+{p)[W{p,p')ij{p')]A^{p) 
3 Ar(p) [W{p,p') ^{p')] ^tip) (17) 



The equivalence of eq.(p^ and eq.(j^) can be shown by 
applying the projectors IT] to both equations from both 
sides. From eq.(|lq) one obtains e.g. 

Atip)m^tip) = o 

Ar(p)V(p)A2-(p) = (18) 

as in eg. (PI) due to the relation Af{p)j'^ = ■j'-* Af{—p). 
Note that eq. dlS) can be also written in the concise form 



LOi UJ2 



(19) 
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B. Normalization condition and scalar product 

The normalization for general BS-amplitudcs has been 
given for bound states with conserved quantum numbers 
by NishiUma Q and Mandelstam We follow Cut- 
cosky who treated the more general case where no 
current has to be conserved. As this has already been 
treated within textbooks (e.g. [Q), we will only give 
the result for the normalization. Let the bound state 
be normalized as (P \P') = {2n f 2P° S^{P - P'). Then 
the normalization condition in graphical representation 
is given by Fig.^ (compare e.g. ||2^). Contracting with 
the momentum of the bound state, this reads explicitly: 



(27r)4 (27r)4 



tr 



dPt' 



iK{P,p,p') xp(p') 



The summation over color indices is suppressed here. / 
denotes the product of the inverse quark-propagators: 



/(P,p,p')c.o',/3;3'=<5W(p-p') i2ny 



(21) 



Note that the vectorial condition of Fig. (^) and the scalar 
normalization (po| ) are in fact equivalent which follows 
from: (a) the formal covariance of the equation and (b) 
the fact that in the rest frame the time component of 
Fig.(^) gives eq.([20[) and the space components vanish, 
as the derivative d/{dP^)[I + K] is proportional to ^p'^ 
or 7* so that the integrals or the trace on the rhs. of 
eq.(|2^) give zero. For an interaction- kernel, which is in- 
stantaneous in the rest frame, i.e. of the type of eq.(|^), 
we have: 



(22) 



so that the contributions of the interaction kernel to 
the normalization vanish. At this point we would like 
to mention that the BS-equation and the normalization 
condition for an instantaneous interaction may be formu- 
lated covariantly, so that the corresponding amplitudes 
X are correctly normalized in any frame. The explicit 
normalization for the corresponding Salpeter amplitudes 
$ will be performed in the rest frame. First we 

define the vertex functions: 



Tp{p) := [S![p^)]-'XP{P) [S![-P2)]-' 
Vp{p) := [S[i-p2)]-'xpip) [S[{pi)]-' 



(23) 



With the BS equation (||) we then obtain the following 
important result: 



(p) P=(M.o, = = 



d^p 

(2^ 



[Vip,p')<i>ip')] (24) 



i.e. the vertex-function depends only on the relative 
three-momentum p. This formula allows the reconstruc- 
tion of the vertex function F and therefore of the full 
BS-amplitude x h'om the Salpeter amplitude Insert- 
ing eq.(p3|) into the normalization condition ( pO| ) the de- 
pendence on p^ is completely determined by the quark- 
propagators, so that the p'^-integration may be performed 
analytically. 

We use the general relation between the BS-amplitude 
X and it's adjoint x for spin-l/2-fermions (see e.g. 
for the scalar case): 



1 



X{p) = -Tr- / dq'' 



(20) xip) 



2711 

1 

2TTi 



f{q°.P, 



9{q°,p) 



q^ + ie p^ — q'^ — ie 



(25) 



7 



( Pil^Pl ^ 9Hl":P) ^ 

^ \pO-q" + ie pO-qO-^' ' ' 



with matrix valued functions / and g. From this we 
derive the following relations in the special case of an 
instantaneous interaction: 

f(p) = -7ort(p)7o $(p)= -/o^Hp)1o (26) 

This leads to the normalization condition for the 
Salpeter-amplitudes in the rest frame: 



(2 



Hp)KipMp)^ti-p)}- (2^)' 



(27) 



2M 



It also may be expressed in terms of the 2x2 amplitudes 
$++ and $ defined in eq.(|l^) as: 



d^p 



2ljJ\L02 



(27r)3 LiJ\m2 ~V uj2m_ 



trl ($++(p))t$++(p) (28) 



- {^—{pyy $""(p) j = (27r)^2M 

This form also shows the connection to the nonrela- 
tivistic norm: in the NR limit <f>++(p) i(T2 becomes the 
usual Schrodinger wave function and $ {p) goes to 
zero as p Vto^ • $++(p). Furthe rmore the weight func- 
tion becomes equal to unity so that we obtain the usual 
Schrodinger normalization. For deeply bound states how- 
ever we have appreciable deviations from this norm, as 
the lower amplitude $ is o f the same order as the up- 
per component $++ (see Sec.IIC). 

Eq. (g^ motivates the definition of a scalar product for 
amplitudes ipi — $17° and ip2 — *i'27° as 



(^1 IV'2) = / tr (v4A+V^2A2--^tAr^2A+) 



tr 



■Ipl {—1p2-1p2 — 
UJl UJ2 



(29) 



with all quantities depending on p and the notation J = 
J d^p/ (27r)^. Note that this scalar product is not positive 
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definite. The normalization condition (|2j) for solutions 
of the Salpeter equation is then given as 



(30) 



The following discussion will be restricted to amplitudes 
satisfying eq.(|l9|). In that case one has 



(V-I W^'^) (31) 



tr 



where the prime indicates the dependence of 'ip2 on p' . 
If one considers kernels that fulfill J j' tr (-01 — 
/ /'tr(-02 Wip'i)* , which is valid for a wide class of in- 
teractions (e.g. for W^}' = / ((p-p')^) Ti V(p')r2 
with hermitian matrices Ti and a scalar function /), the 
Salpeter hamiltonian Ti. is selfadjoint with respect to the 
scalar product given in eq.(p9h, i.e. 



(V'l |HV'2> = (WV^i 1^2) (32) 
which has two important consequences, namely 

• bound state masses M are real for eigen functions 
tp with nonzero norm (■0 1 1/)) 

• amplitudes ■01 and ip2 corresponding to differ- 
ent eigenvalues Mi 7^ are orthogonal, i.e. 

(^1 I ^2) = 

The first point can be seen easily from M (0 | -0) = 
(0 \n ip) = {Hi^ \^) = M* (0 I 0) and the second point 
follows from M2 (V'l | V'2) = {i^i \H 02) = (H0i | 02) = 

Ml* (01 102). 

C. Structure of the Solutions 

The Salpeter equation exhibits some further general 
structures connecting solutions with positive and nega- 
tive eigenvalues. For the case of equal quark masses J.F. 
Lagae has shown that for kernels satisfying [Wtpy = 
Wip'' the eigenvalues will come in pairs of opposite sign, 
the corresponding eigen functions having normalizations 
with opposite sign. In the following we will extend this 
result to the general case of unequal quark masses and 
compare the block structure of the conjugate solutions. 
We further show that nondegenerate bound states with 
mass M — have zero norm (0 | 0) = 0. The discussion 
of physical acceptable solutions is postponed to the end 
of this section. 

To deal with the unequal mass case we first investigate 
the structure of the BS-equation under charge conjuga- 
tion. The details are shown in the appendix with the 
result that solutions of 



(^/2/lV'/2/l)(p) = M^f^f^{p) 



are related through 



(34) 



with Sc given in eq.(D5) in the appendix. The indices of 
Ti denote the flavor dependence of Hi and in eq. (^7|) . 
For simplicity only the case without flavor mixing is con- 
sidered, the generalization being straightforward. 

In the following discussion we will assume that the 
BS-kcrncl satisfies [Wf^^f^^'^^j^)'^ — (0/i/2 )^ which 

is fulfilled e.g. for kernels of the form 1^/1/2^/1/2 ~ 
/{{p — p')'^) r0ji^2(p')r with a hermitian matrix F. 
The hermitian conjugate of eq.(lq) thus leads to 



/2 



(35) 



Renaming /i /2 and comparing this equation to 
eq.(B3) we thus have shown that 



• solutions of the Salpeter equation come in pairs 
(V'/i/^, M) and (V'}^^^, -M*) where eq.(||) con- 
nects the two solutions. 

Consider the normalization of 0/^/2 and '^p\.^f-^'- With 
eq.(|3|) and the relation Sc f^Hp) Sc = *A?(-p) one 

finds (V'/2/i IV'/2/l)/,/i = (^/l/2 IV'/l/2)/,/2 

indices of the scalar product determine the flavors of 
and A J in eq. (p9|) . On the other hand cyclic permu- 
tation under the trace shows that (0^/J V'^/i)/2/i = 
-((^/2/i)^l (^/2/i)^)/iA that we have 

( V'/i /2 1 /2 > /i /2 = - ( V'/, I V-/^ ) /i /2 (36) 
and we flnd that 

• the states with eigenvalues M and — M* have op- 
posite norm. 

Let us now compare the block matrix structure of the 
two conjugated solutions. With the angular momentum 
decomposition eq.(p5|) (see appendix) 



$++(p) 
{P) 



E 

L S 

E 

L S 



(37) 



with the 2x2-matrices 9500 = fiq = <Jq/\/2 it is 

straightforward to show that 



JMj 



/2/1 

E(-i)'^' 



= (-1) 



-,J~Mj 



T^LsiP)] . . [YL{^p)®^stM, (38) 



LS 



hf 



and similarly for $ . It is shown in the appendix that 



[<sb)]/2A = (-l)^+^[^rsb)]/x/2- Since L and S 
(33) are integer the phase vanishes and we obtain the result 
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/ \ JM, 



= (-1) 



J -Mj 



According to eqs.(|l3D,(|14|) we write 



(39) 



(40) 



where the indices of $ indicate the flavor dependence 
of Ci in eq.(|2|). The hermitian conjugate of eq.(|l^) 
with /i and /2 interchanged gives Ajr tp'j:^j^ A^^ = and 
At iIA f At = so that we can write 

II ^ J2jl J2 



(41) 



with some amphtudes C^^, £, that can be determined 
by comparing eq.(^) to 

eq.(g|) with the resuh (C/,/J"'*^-^ = "[(^/s/i )'""^'']^ = 
(— 1)^^'^^*^' (^/i/,)"' ^^'^ ^'^d the same expression with 

++ and interchanged. According to eq.(^l|) we thus 

can write 



(42) 



f-Mj 



We thus have the result that 

• exchanging the functions $++ and <i> in $ = 
$ ($++, $^^) turns an amphtude with eigenvalue 
M into an amplitude with eigenvalue —M*. 

With the relations obtained above it is easy to inves- 
tigate the eigenvalue M — which is assumed to be not 
degenerate apart from the trivial degeneracy in the an- 
gular momentum projection Mj. From eqs.(p3|),(|35|) we 
have 



(43) 



which through eq. (f42| ) implies {''Pf^'fi'V = '^''J-'fih^' '^ith 
|A| = 1. Then eq.(pq) gives 



JM, 

/l/2 



/l/2 



i^iJ'' Hit'') ^^^^ (44) 



Since the scalar product is invariant under rotations one 
can substitute —Mj by AIj and obtains 



JMj 

/l/2 



JMj 

/l/2 



/l/2 







(45) 



So we find that 



• nondegenerate eigen functions with eigenvalue 
M = have zero norm. 



From eq.(|4|) it is clear that (V'/J/)^ = ^^jjjt^' is 
equivalent to setting 



(46) 



This equation illustrates a common aspect of the Salpeter 
equation: in the nonrelativistic limit with M « mi + 
7712 the large component $++ dominates over the small 
component but if one goes to deeply bound states 

(e.g. by increasing the coupling constant of an attractive 
interaction) the two components become more and more 
equal until finally $++ = ±<I> is achieved for AI = 0. 

From the discussion above it has become clear that 
we have to identify the physically acceptable solutions. 
There are two criteria making sure that a solution is ac- 
ceptable: 

• The norm of the solution has to be nonzero which 
automatically implies that M is real. 

• The eigenvalue M and the norm have to be posi- 
tive in order to fulfill the normalization condition 
(V' I V') = (27r)2 2M. 

We would like to mention that the typical doubling of the 
physical eigenvalues is well known from the RPA equa- 
tions in nonrelativistic many particle theory [^ . For 
a hamiltonian H with spectrum En the RPA equations 
have solutions En, —En- This doubling can be traced 
back to the appearance of the time ordering operator T 
in the definition of the particle-hole propagator. There- 
fore neglecting negative mass eigenvalues is consistent 
with the RPA structure of the Salpeter equation (com- 



pare chap. Ill 



The role of the solutions with M = and {^p = 
is not clear. On one hand there is a priori no contra- 
diction with the normalization condition. On the other 
hand the Salpeter equation has been obtained in the rest 
frame of the bound state, i.e. one first performs the limit 
P ^ and then investigates the case M ^ 0. How- 
ever, the correct procedure for massless bound states is 
first to perform the limit M ^ in the BS-equation. 
In the resulting equation one then can study the limit 
P ^ 0. It cannot generally be expected that exchanging 
the limits for P and M leads to equivalent results (com- 
pare 01 for a more detailed discussion of this problem). 
Furthermore the definition of the instantaneous interac- 
tion kernel eq.(|6|) becomes dubious since p± is not well 
defined for P ^ 0. We therefore prefer to require M > 
and ("0 I V') > for physically acceptable solutions. 



III. NUMERICAL TREATMENT 



From the definition of $ in eq.(|13[) it is easy to derive 
a basis expansion for tp — $7°. Let 



M, 



(47) 
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be a complete set of 2 x 2 basis functions with real radial 
functions RmLiip)- The basis functions are chosen or- 
thonormal with respect to the usual scalar product given 
by 

d3p 



which is of the same form as the well-known RPA equa- 
tions in nuclear physics jl^]. Let (a^+\ a(_2) be an eigen- 
vector with eigenvalue M. Then eq.(54) shows that 
{{a^~^)*, (a^^^)*) is an eigenvector with eigenvalue —M* 
which is just the result of the previous section. Usu- 



{Ei \Ej)= f -A^tr [E+{pjEj{p)] =% (48) ally the Salpeter hamiltonian H has the property that 
J i^'"') the matrix elements Hff and also the eigenvector coefH- 



where the trace just gives the usual scalar product for 
the spin matrices ti(p^j^^(ps'M'^ = Sss'^MsM'^- Note 
that the angular structure of the basis function s ma tches 
the structure of $++ and as given in eq. (|B5|) . We 
choose RmLiip) to be real functions. It is now possible 
to expand 



(49) 
(50) 



which implies TZ^^siP) = af^' RnihAp) ShhiSsSi 

for the radial wave function. Since Ti-^^g (p) are real func- 
tions in most cases of interest, the coefficients a[^^ then 
also have to be real. Now define the 4x4-amplitudes 



(±) 



and also the eigenvector coeffi- 
cients a[^^ are real within the basis given above. In that 
case and since Nff is real, also M must be real. This 
result has already been shown before for eigen vectors 
with nonzero norm. 

Furthermore we see that if M = is an eigenvalue we 
expect from eq.(]5^) that the eigenvector fulfills = 
±a'~^ and we reobtain the result that this solution has 
zero norm. 

In a numerical treatment only a finite basis i < imax ~ 
10 can be taken into account. Then eq.(^J) becomes a 
finite matrix equation that can be solved with standard 
numerical methods. One thus obtains an approximate 
eigenvalue Mp and an approximate eigen function ipp to 
the Salpeter equation that exactly fulfill the relation 



(55) 



=(+) 

S-) _ 



$(£;., 0)7° 



= <i>(0,i?,)7° 



(51) 



Note that these functions are not orthogonal with respect 
to the scalar product given in eq.([29|). Since <l is bilinear 
we nevertheless can expand ?/; as 



The index (3 indicates that the basis states Ei and thus 
ipp depend on a variational parameter j3 of dimension 
MeV~^ that sets the absolute scale for the momentum 
dependence via E^{p) = P^''^ Ef^\pP). 

If solutions of nonzero norm are considered the 
Salpeter equation is equivalent to the variational prob- 
lem i 



(52) 



SM[^P] = S 



= 



(56) 



_ _ where the variation 6 is taken over all functions tp with 

so that the constraint i^A^ = ijjA^ = is auto- nonzero norm that fulfiU A+^A+ = A'^pA' = 0. Ac- 



matically fulfilled. The Salpeter equation Htp = Mtp can 
now be written as the matrix equation 



H- 



H- 



= M 



N- 



N+- 
N— 



(53) 



with Hlf = {ef'\rLeY') and iV//' = (e,^"'^|e^"'^). From 
the definition of the scalar product one easily sees that 
N++ = -N — and N+- = N-+ = 0. Furthermore 
we find from eqs.(^),(|l]) that e|^' and e\ ^ are con- 

nected by [{e^Xh^^ = (-l)^"^"^^ (e!-));?jf ^ so 
that we can use eqs.(|35|),(p6|) as well as the invariance 
of the scalar product under fiavor exchange and under 
the replacement —Mj Mj to obtain H^j~ — {H^^~^)* 



cording to eq. (|55|) we make the variational ansatz ijj = ipp 
implying M[tpp] = M(j and look for stationary points of 
Mp as a function of (3 (which also fixes (3 for each meson) . 

The calculation of the matrix elements within an ex- 
plicit model will be postponed to our second paper. At 
this point we would only like to make a few technical 
comments. The matrix elements of the interaction ker- 
nel can be efficiently calculated by inserting two complete 
sets of basis functions written schematically as 



{^\h{p)V{r)h{p')\j) 

= E ^9\V{r)\h) {h\h{p')\j) 



(57) 



and H- 



The matrix representation of the 



Salpeter equation thus takes the form 



H++ H+- 
{H+-)* {H++) 

N+- 




= M 







a(-) 



(54) 



so that V{r) can be parameterized in coordinate space. A 
suitable choice for the basis functions is given by the func- 
tions RnL{y) = N^L Ll^-^Hv) e-yl"^ with y = pf3 and 
L'^^'^iy) being a Laguerre polynomial. We found that 
about ten basis states are sufficient to solve the Salpeter 
equation with rather high accuracy. The choice of 3- 
dimensional harmonic oscillator functions is less favored 
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since their asymptotic behavior ^ for y — > oo 

turns out to be not appropriate for our quark model (es- 
pecially for deeply bound states like the Pion). 



IV. DECAY OBSERVABLES 

Apart from describing the mass spectrum of mesons, 
any realistic model must also be able to describe mesonic 
transitions and decays. The important question arises 
whether a good description of the extremely deep bound 
states as pion or kaon can be combined with a reasonable 
description of confinement. The BS formalism offers a 
natural framework, as the role of the lower component of 
the wave function turns out to be crucial for the correct 
normalization and calculation of the decays. This can 
be seen most clearly in the following formulas for the 
leptonic decays and the weak decay constants. 

We furthermore show how to reconstruct the full BS 
amplitude, which may be used as a starting point for the 
calculation of electromagnetic or hadronic transitions. 



A. Leptonic decay width and weak decay constant 

The transition or decay of bound states are calculated 
from BS amplitudes using the formalism given by Man- 
dclstam . We will merely sketch it by considering first 
the leptonic decays of vector mesons. The corresponding 
Feynman diagram is given by Fig.^ where Kiep denotes 
the kernel which is irreducible with respect to the incom- 
ing qq-pair and the outgoing /"'"Z~-pair. If we consider 
only graphs of leading order in the electromagnetic cou- 
pling constant we obtain the approximation on the right 
hand side. 

The exceptional role of these decays (together with the 
weak decay constants) is that if the BS kernel would be 
exact also the decay amplitudes would be correct to any 
order in the strong interaction. The hadronic part of the 
transition matrix clement is given by p9[ |: 



(0 b-^(o)| P1-) 



■tri^^^p)) (58) 



The instantaneous approximation thus simplyfies this 
calculation, as one can express the transition in terms of 
the Salpeter amplitude It may be decomposed using 
the transformation properties under rotation and par- 
ity, as is shown in appendix A (compare eq.( |E2| )). For 
= 1~ mesons it reads (with p := |p |): 

<^>'"'{p ) ^ ^ToiP) YiMi^^p) + ^%{p) 70 yiAf (^^p) (59) 
+'^11{P) 7570 [Yi{np) X 7]^^, + $iO(p) 75 [Y,{np) x 7]^^, 

-H$[!°(p) [Yoin^) X 7]^,, + $[!}(p)7o [Yoinp) x 7]^^, 



The radial amplitudes are determined completely 



by the functions $++ and $ 
latter have the form: 



For a 1 meson the 



±± 

IM 



+ ni'^\p) [i2(r!p) X ^1], 



(60) 



The integration and trace in eq.(pq) for the current pick 
up only the s-wave amplitude The usual spin sum- 
mation and averaging leads to the decay width: 



r(r 



24 



M3 



l+r) = (61) 

0(<'(p)-KiT'(P)f 



where a = 1/137 is the electromagnetic coupling con- 
stant and Cq is the fraction of the quark charge com- 
pared to the electron. Apart from the contribution of 
the lower component <I> to the normalization eg. (p8|) , 
its importance comes out here very clearly. For M we 
use the experimental meson mass to obtain the correct 
phase space. 

The next observables to be considered are the weak 
decay constants /jr and . They are defined by the ma- 
trix element of the axial current pl| ] (with this definition 
j.{exp) ^ -^22 MeV): 



(62) 



Again for an instantaneous interaction this can be eval- 
uated from the Salpeter amplitude 



(63) 



and in terms of the 

73 f (fp 



T^oa ^ Yqo (poo- 




ip) - 7^, 



(-) 

00 



(64) 



B. The decay tv° , 



27 



These decays provide another test for the description 
of the low lying pseudoscalar mesons. To our knowl- 
edge they have not been calculated in the framework of 
the full Salpeter equation (for a slightly more restricted 
ansatz see Mitra et al. [p2|). The basic idea here is to 
reconstruct the vertex function F from the Salpeter am- 
plitude $ by means of the BS-eq. ( p4[ ) itself, which gives 
the full fourdimensional structure of the BS amplitude x- 
This has to be taken into account for a correct descrip- 
tion of decays. The corresponding Feynman diagrams for 
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the neutral pseudoscalar decay in lowest order of the in- 
teraction are given in FigJJ. Of course this is not correct 
to any order in the strong interaction like e.g. in the case 
of the pion decay constant, as we obviously neglect the 
strong interaction of the intermediate quark. Therefore 
we expect these calculations to be less accurate. 

The T-matrix element in lowest order for a meson with 
mass M decaying into two photons with momenta fci , k2 
is given by 



T 



zV3(ieg) 



trl {P/2+p)T{p ) {^P/2+p) 



^1 5^(-P/2+p-fci)^2 



5^(-P/2+p-fc2) 



(65) 



Therefore we have to reconstruct the vertex function 



r from the Salpeter amplitude $ by means of eq.(24). 
Both of th em have an angular decomposition according 
to eq.(E2). The four scalar amplitudes depending only 
on \p \ are expanded in a basis of radial wave functions. 
The quark propagator can be written as 



A+{p) 



A-ip) 



Po - u + ie po +Ll! ~ ie 



(66) 



which is suitable to perform the g° integration using the 
residue theorem. Then we substitute p —p in the sec- 
ond term of eq.(|65|) and use T{—p) = — 7°r(p)7° for 
pseudoscalar mesons which follows from eg. (p5[) . We 
use {7*^,^^}+ — for transversal photons and define the 
quantities 

e± (fc -p):= ^2 A± {k - p )7Vi - ^1 A± (fc - p )772 (67) 

with ujk := {{k — p Y + m?)^/'^. Furthermore the product 
of the vertex function with the projection operators 
is denoted as: 



r+-(p) := A+(p)7"r(p)7°A-(-p) 



(68) 



and r++ etc. analogously. In concrete calculations these 
quantities are evaluated using a matrix formalism, which 
represents the multiplication of T with a Dirac matrix 
by a transformation amongst the scalar amplitudes (see 
(E4)). After some calculation we obtain: 



(Pp 
(2^ 



E 



(69) 



(A/ /2±LJk± uj) (-M/2 ± Wfc ± w) 



rT± Q^{k^p) 



{M ±2uj){Ml2±uik±oj) 



{M ±2uj){M/2±ujk±i 



where the summation ± runs over the upper and lower 
sign. For Q we obtain 



e±(fc - p) = i{ei X £2) • S*(P - fc) 



(70) 



with 



E±(p-fc):=75 j^^± -l(_/(p-fc) + TO707)^ (71) 

To avoid the zeros in the denominator for Ad > 2m in 
eq.(|6|) we use the BS equation in the form ^ to obtain 
finally (with fp = t+ -f.-): 



T = V3 (£1 X £2) ■ S{k) with 
(fip 1 



S{k):^ 



(72) 
(73) 



(27r)3 M2/4- (c^fe +cj)2 
tr|r++(p) t-{p- k) - V—{p) t+{p- k) 
M 



M + 2oj 



r-+(p) Y}>{p-k) 



i M 



As for transverse photons we have £1 _L fc _L £2, for 
k — kcz the only nonvanishing component of £1 x £2 
is in z-direction and we only need to calculate the z- 
component of S. The standard formula for the decay 
rate of a particle with mass M then yields: 



r(7r0,ry->27) 



3 «2g4 



167r M 



^|5.(fce-)|^ 



(74) 



V. CONCLUSION 

We have investigated the structure of the instanta- 
neous BS-equation (Salpeter equation) for the general 
case of unequal quark masses. Furthermore we have de- 
veloped a numerical scheme to solve the Salpeter equa- 
tion which enables the calculation of mass spectra and 
Salpeter amplitudes. In order to test various models 
beyond the mere reproduction of the mass spectra, we 
have further given explicit formulas for the computation 
of weak meson decay constants (/,r, Jk etc.), the decay 
widths into two photons and into an electron-positron 
pair. 

Because of the relativistic kinematics, the correct rela- 
tivistic normalization of the amplitudes and the dynami- 
cal treatment of the lower component $ we expect the 
Salpeter equation to provide a framework for quark mod- 
els that is superior to other treatments like the reduced 
Salpeter equation or the nonrelativistic quark model. We 
will investigate an explicit quark model based on the 
Salpeter equation in a subsequent paper |l|. 
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APPENDIX A: SPECIAL LORENTZ 
TRANSFORMATIONS 

Let A be a special Lorentz transformation and g be the 
corresponding element of the covering group SL{2,C), 
given by 



where 



g (t{x) .gt = (j{Kx) 



(Al) 



with a{x) — x^cTf^ and (cr^) = (!,(?). The transformation 
matrix Sg for Dirac spinors in the Weyl representation is 
then given by 



Sn 



9 

(g^y 



With the transformation matrix 
B 



1 



1 1 

1 -1 



(A2) 



(A3) 



we can use — B B^^ to transform Sg into the 
standard Dirac basis employed in this work. The result 
reads 



s -1 



9 + (9^)-' 9- (9^)-' 
9~{9^)-' 9+ [9^)-' 



(A4) 



Note that the relation A'^^, •y'^ — Sg^ Sg is automat- 
ically fulfilled. For a boost A with A (M, 0) = P and 
~ one has explicitely 



M 



2{M + PO) 



(A5) 



and for a 3-dimensional rotation A = P one has g = u (z 
SU{2) and therefore 



'S'ti 



u 
u 



(A6) 



in the Weyl basis as well as in the standard Dirac basis. 

The transformation of a Dirac field operator '^{x) is 
given by 



Ug^{x)U-' =S;'^iAx) 

Ug^ix)U-' =^iAx)Sg 



(A7) 
(A8) 



with Ug being a unitary operator. Consider a 
state with momentum P and P^ = M^, total an- 
gular momentum J and 3-component Mj normal- 
ized as {P',J',M'j\P,J,Mj) = (27r)3 2pO(53(p' _ 
P) 5j' J Sm'jM J ■ Then the action of Ug on this state is 
given by 

Ug |P, J, Mj) = I ^'^'j) ^M', M, («) (A9) 



1/2 

is a Wigner D-function with D^^, ^^^^ (u) ^ um'^ Mj and 
u = [cr(AP/A/)]i/2 g [cr(p/Af )]i/2 is the Wi gncr rotation. 

From the definition of the BS-amplitudc in Sec. II A one 
can see by inserting the unity operator 1 = U~^ Ug that 
the BS-amplitude transforms as 

ip) - E ^ip'' (^P) ^9 Dir, M, (") (AlO) 

The BS-equation is compatible with this transformation 
law for covariant kernels. 



APPENDIX B: ANGULAR DECOMPOSITION OF 
THE 2x2- AMPLITUDES 

Let A = P be a 3-dimensional rotation and q ^ u be 
the corresponding matrix g SU{2). With eq.( Al(]| ) and 
the block matrix structure of $ given in eq.(|lCl|) as well 
as the relation (—iao) u^^ iao = *u we find 



[^++mJM,icj2 (Bl) 

and the same for the other amplitudes where we 
use the tensor notation {[u x u] {^~^^ ia2)}mim2 = 
J2ni[7n'^ Ur,n7n[ Um2Tn'2 {^'^^ i<^2),n[m'2- Therefore One Can 
identify the indices mi, m2 of the matrix elements 
[$++ i(T2]^^„j^ with the 3-components of the spin of 
the quark and the antiquark. Let {xsMs)mim2 — 
{l/2mil/2m2\SMs) be the spin matrix of the two 
quarks coupled to the total spin S. Define ^psqi<^2 
XSq, i-e. 



or explicitely 



1 



1 



1 

71' 



1 

71 (o 



1 
1 



-1 



",511 



-1 




1 



Then eq.(Bl) implies that we can decompose 
as 



L S 

E 

L S 



(B2) 

(B3) 
(B4) 

(B5) 



with the spin S and the orbital angular momentum L 
coupled to J. We assume that the BS-kernel allows 
choosing Tl^^g{p) to be a real function. The sum goes 
over all values L, S that are compatible with parity and 
charge parity of the bound state, see below. 
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APPENDIX C: PARITY TRANSFORMATION OF 
THE BS-AMPLITUDE 

The Dirac field operator ^(a;) and the bound state 
with parity number ttp = ±1 transform under parity 
transformation P as 

Up Up^ = 7" *(Px) (CI) 
C/p*(x)C/p^ = *(Pa;)7° (C2) 
Up \K, J, Mj, Tip) = TTP \K, J, Mj, ttp) (C3) 

where K is the momentum of the bound state and P x — 
{x^,—x). For the BS-ampHtude this imphes 

XK(p) = 7rp7"xp^(Pp)70 (C4) 

For the block matrix structure one finds 

<i>-(p) ci>-+(p) j -''^ 

(C5) 

To be compatible with the angular decomposition 
eq.(B5) we find the well-known condition ttp — (—1)^+^. 



APPENDIX D: CHARGE CONJUGATION 

Let 1/1/2, P) be a qq-bound state with flavors /i and 
/2 and momentum P. The charge conjugation then acts 
like 

C/c*a(x)C/c^ -^(5c)a/3**^(x) (Dl) 

13 

Uc'K{x)Uc' ^Y.^ScU^0{x) (D2) 

Uc\hhP) = \f2h,P) (D3) 
and we find for the BS-amplitude 

xifM = -iSclo) 'xmA~p) (^c7°) (D4) 



with the matrix 5(77*^ given in the standard basis as 
Making the choice 771 = 772 = 1/2 implies 



*7J. <f 7i y I *7r/. (-^^ "i^r/. (-^^ ; 

(D6) 

with $++ = —ia2 (*<1'^^) jo'2- With the angular decom- 
position eq. (|B5| ) and ifis q ~ fSq one finds 



/1/2 



/2/1 



(D7) 



For an eigenstate \P, ttq) of the charge conjugation we 
have /i = /2 and Uc \P, t^c) = t^c \Pi t^c) which implies 
the well-known condition ttc = (—1)^+'^. 

To investigate the compatibility of eg. (p4|) with 
the BS-equation we use the relation Sjlj>i) = 
— (5*07°) '^Sj,{—pi){ScJ^) for the fermion propagator. 
Let Xfif2{p) be a solution of the BS-equation with /i 
being the quark flavor and /2 being the antiquark fla- 
vor. It is straightforward to show that Xf2fi{p) given 
in eq.(D4) is a solution of the BS-equation with inter- 
changed flavors. For the Salpeter equation this implies 
that the solution "0/1/2 ^-^d the solution of the equation 
with flavors interchanged 0/2/1 ^-re connected by 



(D8) 



APPENDIX E: ANGULAR DECOMPOSITION OF 
THE 4x4 AMPLITUDE 

From the angular decomposition of the 2x2 blocks and 
the property under parity transformation we also obtain 
the structure of the 4x4 amplitudes. 

With the spherical components cr„,. of the Pauli matri- 
ces we define the two tensors 
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[0] 



1 
1 



inis 











(El) 



which are again Dirac matrices. We obtain the following 
representation for the BS amplitudes in the rest frame 



(E2) 



X 7^' 



Mj 



The sum over i covers all the values of Si,ai,gi e 
{0, 1} and li G {J, J ± 1} that are compatible with 
(^_iyi+h+ai _ where up is the parity of the meson. 
In general we have eight scalar amplitudes Xi{P^ ^ \p\) (ex- 
cept for mesons with spin 0, where there are only four). 
The same structure holds for the Salpeter amplitudes 
the vertex functions F and any amplitude that is built 
out of these by multiplication with e.g. 7°, 7!]^ or A^(p ). 
It is therefore convenient to introduce a matrix formalism 
for such multiplications which is useful for the calculation 
of decays or tra nsitions and simplifies the Dirac algebra. 
From eq.(E3) we may write the amplitude in the form 



(E3) 



Hp) = E '^'^(IpI) 



i=l 
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so that we can express the muhiphcation from the left or 
right hand side with a matrix 7 e {7*^, -fp, ...} by a trans- 
formation Dl{j) or -Dij(7) in the 8 dimensional space of 
the radial amplitudes: 



(E4) 



i,fe=l 



and analogously with Djt for multiplication from the 

right side. We also use the following formula to calcu- 
late the trace in Dirac space and the angular integration: 

Jdflp tr { rliQp) } = 4 MW\) (E5) 





FIG. 1. Graphical representation of the BS equation 



+ 



iK 




= 2iP, 



FIG. 2. The normalization condition 
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FIG. 3. Leptonic decays in the Mandelstam formalism up 
to lowest order in the electromagnetic coupling constant 





FIG. 4. The decay 7r°,?? ^ 27 
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